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The implications of the spin exciton mechanism are exposed in the context of a Spin Density 
Wave (SDW) instability occurring inside the superconducting phase of a layered heavy electron 
compound. In this model a magnetic field serves as a tuning parameter bringing the system to the 
point where the transverse spin correlations are enhanced due to d I 2_ H 2-wave superconductivity 
and induces an instability to a phase with coexisting superconductivity and SDW order. The 
model considers electrons in a crystal with antiferromagnetic interactions and provides restrictions 
on the Fermi surface characteristics and on the ordering wave-vector (that can be commensurate 
or incommensurate close to commensuration). The applications of the model are addressed to the 
low-temperature/high-magnetic-field phase of CeCoIns [M. Kenzelmann et al, Science 321, 1652 
(2008), Phys. Rev. Lett. 104, 127001 (2010)P 
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A significant experimental finding was achieved by 
Kenzelmann and collaborators^ while using Elastic Neu- 
tron Scattering to probe magnetism in the supercon- 
ducting phase of CeCoIns at low temperature and un- 
der a strong magnetic field applied in the basal plane 
of its tetragonal crystal structure. The upper criti- 
cal field in CeCoIn 5 is mostly determined by paramag- 
netic limiting (H c2 (T = 0) ~ 11.7 T) and due to this 
the phase transition to the superconducting state be- 
low T = 0.4T C (T c = 2.3 K) 2 is of the first ordeM In 
papei^ the authors observed a magnetic ordering with 
wave- vector Q = (q, q, 1/2) where q ~ 0.45 is an incom- 
mensurate wave vector (here we use Reciprocal Lattice 
Units in which 2n/a = 1, o being the in-plane lattice 
constant), which was independent of the field magnitude. 
The small value of the magnetic moment on Cerium sites 
to = 0.15/is (pb is the Bohr magneton) oriented along 
the c-tetragonal axis indicated that the ordering is of itin- 
erant origin. It was remarkable there that incommensu- 
rate SDW was confined inside the superconducting phase, 
meaning that here superconductivity is an essential ingre- 
dient for SDW to develop. The existence of the magnetic 
order was first detected by the technique of NMRpB and 
its precise field-dependence later determined^. 

Theoretically, these coexisting orders were discussed 
within mod els that can be divided into three classes. In 
the first one^H2l theories rely on coupling between SDW, 
superconductivity and superconductivity with Cooper 
pairs having non-zero center of mass momentum (Fuldc- 
Fcrrel-Larkin-Ovchinnikov (FFLO) phase, Pair Density 
Wave (PDW) , or 7r-triplet superconductivity) , which sta- 
bilizes SDW and superconductivity at high-field and low- 
temperature. Evidence for the occurence of a staggered 
superconducting state in CeCoIns is however still to be 
revealed. A second point of view^ highlighted the role 
played by the vortex lattice which can increase the den- 
sity of states in the nodal direction o f the gap and trigger 
a magnetic instability. In the thirclUl^ ft was empha- 
sized the importance of Pauli limiting in d-wave super- 
conductors for stabilizing SDW order in the case where 



the ordering wave-vector is the nesting wave-vector Qn 
(in the sense that it joins two points of the normal- 
state Fermi line where the Fermi velocities are paral- 
lel). There was found that low-temperature/high- field 
superconducting phase in CeCoIns is a coexisting phase 
of FFLO and incommensurate SDW orders.^ Also, there 
was pointed oulP^ the enhancement of nesting by super- 
conductivity in the gap nodal direction. 

Another important observation was made on CeCoIns 
thanks to Inelastic Neutron Scattering (INS). Stock and 
collaborators^ measured a spin resonance that was sharp 
in energy (ui = 0.60 ± 0.03 meV) and having wave vec- 
tor distribution centered on Q = (1/2,1/2,1/2) with a 
width of ~ 0.15. Thereafter Panarin and collaborators^ 
studied the evolution of the resonance in a magnetic field 
applied in the [1,-1,0] direction. They observed the same 
resonance with a decrease in its energy and a broadening 
in its lineshape as the field increased. They were able to 
measure it up to ~ 0.5H C 2 where the signal was lost in 
the incoherent part of the spectrum. 

Theoretically Eremin and collaborators^ 5 -' have at- 
tributed the resonance to the proximity to the threshold 
of the particle-hole excitations continuum which is at en- 
ergy uj c = min(|Ak| + |Ak+ q |). Another scenario related 
to a magnon excitation was proposed by Chubukov and 
Gor'kov^ 1 . 

In this paper we present a new direction in the in- 
terpretation of the occurence of a phase with coexist- 
ing superconduc tivity and SDW in CeCoIns where the 
phenomengP an d 13 * 14 ! are closely connected: the reso- 
nance that exists at uj = 0.6 meV shifts to lower ener- 
gies as a transverse magnetic field is applied and triggers 
a magnetic instability before superconductivity is sup- 
pressed. We consider here a situation with electrons in a 
crystal having antiferromagnetic correlations as was pre- 
viously discusse d 13 * 15 ^. 

We use a model of a two-dimensional system and show 
that it provides a consistent scenario for the presence in 
CeCoIns of a SDW order that is confined inside the su- 
perconducting state without requiring the coupling with 
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another state like FFLO or PDW, or nesting proper- 
ties of the Fermi surface. The applicability of a two- 
dimensional model for CeCoIn 5 was already discussed 1 -' 
because of a lack of strict two dimensionality for the 
compound. Here we argue that the Fermi surface in 
CeCoIns has sufficient 2D character so that we can con- 
sider a model of a 2D metal from which superconduc- 
tivity develops. This a ssump tion is corroborated by ex- 
perimental observation d 17 l 18 l In particular,^ points out 
that the direction (1,1,0) manifests particularly strong 
two-dimensionality. 

We found that the c-axis static susceptibility increases 
under a field directed in the basal plane due to the prox- 
imity to the resonance and becomes larger than in the 
normal state. The wave-vector of the SDW is not con- 
strained by the in-plane orientation of the field (as re- 
marked experimentally 1 ) but by the gapped energy spec- 
trum characteristics. The conditions for stabilization 
of magnetism with respect to the normal state at fi- 
nite magnetic field are the following: (i) The Fermi line 
must contain hotspots with ordering wave-vector close 
to (1/2,1/2) corresponding to antiferromagnetic corre- 
lations), (ii) The two values of the gap on hot-spots 
must be non-zero and of opposite signs. As already 
emphasized 15 for CeCoIns, this requires the symmetry 
of the gap to be d x 2_„2-wave. This gives an upper 
boundary for the value of the field that induces SDW 
IJ-Bsdw < (Ahsi + A hs2 )/2, where fx = 5/xb/2 is the 
electron magnetic moment. 

One way to treat at the same time superconducting 
and antiferromagnetic correlations was put forward by 
Scalapincff^l who dealt with Coulomb's interaction in a 
Random Phase Approximation susceptibility %(q, w) = 
Xo(q,w)/[l - UqXoi^ui)], where Xo(q,^)) is the suscep- 
tibility that includes superconducting correlations and 
U q is a momentum-dependent Hubbard Coulomb re- 
pulsion potential (this was originally introduced for q- 
independent interaction). In this model the conditions 
for a collective excitation (called spin exciton) to occur 
are U q Uexo(<l, w) = 1, and C/ g 3mxo(q, w) < 1. 

Here superconductivity is accounted for in the mean 
field BCS theory for electrons under a magnetic field that 
couples to the electron spin via the Zeeman energy (the 
effect of vortices on the spin-exciton are discussed irP2 i n 
the context of high-T c superconductors). The Gor'kov 
Green's functions write once compacted into the Nambu 
(particle-hole space) notation 

g ^. _ (jgn + CT^B)tq + AfcTl + £ fc T 3 

CT (iCj n + a/iB — Ek){iCbn + a\xB + E k ) 

Here Q n — u> n + sign(u) n ) / (2t) , uj n — irT(n + 1/2) are 
Matsubara frequencies, and r is the electron relaxation 
time. We consider a gap order parameter A k which is 
uniform in space and its phase was set equal to zero, 
Ek = + A\ is the zero-field energy of excitations in 
the superconducting state, tq is the unity matrix and T\ 
and r 3 the Pauli matrices. Throughout we take the spin 
quantization axis the same as the magnetic field direc- 



tion which is fixed to be the z-direction in our spin-space 
frame and belongs to the basal plane of the crystal. 

The normal-state electron energy spectrum in CeCoIn 5 
was calculated^ in a model involving hybridized 
conduction-electron bands and f-electron bands with 
Coulomb's interaction. Here we take the two-dimensional 
spectrum 

ik = 2t[cos(fc a ) + cos(feb)] + it' cos(k a ) cos(fcb) 

+2i"[cos(2fc Q ) + cos(2fc fc )] + e, (1) 

with t' = -0.5i, t" = 0.4i and e = 0M. We con- 
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FIG. 1. Contour plot of the normal-state energy spectrum 
= (the Fermi line) in the first Brillouin zone. The con- 
figurations of the commensurate and nesting wave- vectors are 
drawn. 

sider a singlet d^ ^2 -wave superconducting state with 
gap Ak = (A /2)[cos(fc a ) — cos(fcfc)] corresponding to the 
irreducible representation Bi g of the tetragonal crystal 
point group. In the following all energies are counted 
with respect to the scale given by the gap A at zero 
temperature and zero field, and we set t/Ao = 7. 
The free electron spin susceptibility is defined as 

xS b (q, ivm) 

= -- ^2 Tr [ <7 "<7'^'( k )^n)°'a' ( r^o-(k + q,ia;„ + w m )] , 

k,uj n ,cr,<7' 

where the trace is performed in the particle-hole (Nambu) 
space and a° a , are the Pauli matrices in spin space. The 
susceptibility is calculated by taking account of the self- 
consistency equation for the order parameter, yielding 
slight decrease in the gap magnitude at high field. We 
evaluate the retarded susceptibility along the crystal c- 
axis (x-direction of the frame in spin-space introduced 
above) under a magnetic field applied in the crystal basal 
plane (z-direction) 

xS x (q^)= (2) 
\ E { - /2 ( k ' tilfP" - a ^ B ) - K E k+* + v^B)} 

1 

X -Ek - E k+q - 2cr fiB + uj + i/(2r) 
1 

+ E k - E k+q - 2o fiB -u- i/{2t) . 
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-p 2 (k,q)[l-f(E k 



Eu + E. 



k+q 



ofiB) - f(E k+q - a/iB)} 
1 

2a (iB + uj + i/(2r) 
1 



Eft + Ek+q - 
j2l„2 f/n^\2 



2cr /iB 



t/(2r) 



} 



Here ^ k = J B z d 2 ka 2 /(2n) 2 , the integral being evalu- 
ated over the (first) Brillouin zone fc Q ,fcb G [— 7r/a, 7r/a]. 
The coherence factors are 



J 2 (k,q) 



; i 



£fc£fc+g + A k A k+q 



EkEk 



+1 



(3) 



P 2 (M = J(i-^^). (4) 

2 \ &k&k+q J 

Under the condition of a superconducting gap having 
d a; 2_ 1) 2-wave symmetry, the coherence factor p 2 is close to 
unity at the points of the Fermi line where A k = — A k + q - 
The importance of the s ymmetry of the order param- 
eter was emphasize d 22 * 23 ^ for the occurrence of the res- 
onance. Throughout we considered the temperature 
T = 0.05A , the g-factor g = 2, and the damping pa- 
rameter l/(2r) = 0.02A . 
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FIG. 2. Condition for obtaining an instability from the su- 
perconductivity to coexisting superconductivity and SDW or- 
dering at finite field in the case of incommensurate wave- 
vector qic = 0.45. The grey line represents the normal- 
state static susceptibility along the crystal c-axis xW = 
E kCT [/ K* - °I*B) - f (&+, + ofiB)} x - - 2afiB] / [(& - 
£fc+q — 2cr/iB) 2 -|-l/(2r) 2 ] and the black line the superconduct- 
ing one. 

The field-dependence of the real static susceptibility 
Xq x {w = 0,B) is shown in Fig. |2]) for incommensu- 
rate ordering wave-vector and illustrates the condition 
for the magnetic instability to occur, which is Xq x (SIj w = 
0) = 1/Uq with the additional condition for the exis- 
tence of superconductivity B < H p (the zero tempera- 
ture paramagnetic critical field was computed from the 
superconducting free energy with the band structure in- 
troduced above and was found to be H p = 0.36A //i). 
Due to d x 2_ y 2-wave superconductivity the free suscep- 
tibility increases with field, becomes larger than the 
normal-state susceptibility and reaches a maximum at 



the value B = A^ s //i. Before this point a field can possi- 
bly trigger a SDW instability in contrast to the normal- 
state case where Xon(B) remains essentially constant. 
For arbitrary direction, the free spin susceptibility writes 
Xq V = cos 2 (ip)xo z + sin (y>)xjj x , where ip is the angle be- 
tween the magnetic field and the c-axis. The longitudi- 
nal susceptibility Xq z doesn't carry the field dependence 
in the denominator and therefore for B||c no magnetic 
ordering along the c-axis can be induced. This point is 
consistent with experiment!^ where it was found that sig- 
nal of magnetic ordering disappears as the angle between 
the field and the crystal plane increases. 
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FIG. 3. Above: real (dotted lines) and imaginary (full lines) 
parts of the RPA susceptibility computed at magnetic fields 
fiB/Ao = (brown), fiB*/(2A ) = 0.11465 (purple), and 
fiB* I Ao = 0.2293 (green). We used the potential energy 
value Uqic/Ao — 1/0.0583. See text for discussion. Below: 
evolution of the resonance energy with a magnetic field ap- 
plied parallel to the crystal plane. The value of the conden- 
sation field corresponding to vanishing resonance energy is 
found to be B* = 0.2293A //x. 

The upper part of Fig. (|3| shows the evolution with 
field of the collective excitation as given by the RPA 
dynamic susceptibility at wave-vector Q = Qic — 
(0.45,0.45). The value of the gap at hot-spots is A^ s = 
0.27A . We here discuss several points: (i) The zero 
field collective peak appears at energy uj res ~ 0.384A < 
2A/ ls . (ii) Under a transverse magnetic field a well de- 
fined collective peak is shifted to lower energy due to 
Zeeman splitting of the energy of elementary excitations. 
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The higher energy feature of 3mx II (w) is not present 
since the excitations are strongly damped by the contin- 
uum. This observation is consistent with experiment!^, 
(iii) The field at which the excitation softens to zero en- 
ergy is fiB* = 0.2293Ao < ui res /2, hence showing the 
transition from the finite energy resonance (that corre- 
sponds to the condition f/oSftcxofaj w ) = 1) to the ground- 
state instability at the spin-exciton condensation field. 
The evolution of the resonance energy with applied field 
is represented in the lower part of Fig. ([3| . We emphasize 
here the necessity for the transition from superconduct- 
ing to normal state to be of the first order since a vanish- 
ingly small gap would give negligibly small effect. The 
physics behind the effect described here is reminiscent of 
the excitonic phases close to a semiconductor/semi-metal 
transition^. Here excitons originate from d 2; 2_ y 2-wave 
superconductivity under a time-dependent perturbation 
and magnetic ordering represents the zero-energy con- 
densation of excitons at finite wave- vector under a trans- 
verse magnetic field. 

CeCoIn 5 is a system close to antiferromagnetism and 
for this reason the interaction U q is considered to be 
maximum at the antiferromagnetic wave-vector and as 
a consequence the collective excitation dispersion is 
predicted^ to be centered on the commensurate wave- 
vector. The dispersion shows^ a downward shape and an 
incommensurate ground state is expected. This situation 
is however sensitive to the precise band structure, which 
is knowrP^ to be quite complicated and to consist of sev- 
eral bands. Experimental determination in CeCoIns of 



the excitation dispersion would represent progress on this 
issue. 

To summarize, we have presented a new approach for 
understanding magnetism that is tied with superconduc- 
tivity in CeCoIns. The condensation of the induced 
by d x 2_ y 2-wave superconductivity spin-exciton driven by 
the Zeeman splitting of the energy of elementary excita- 
tions generates a state with coexisting magnetism and 
superconductivity. There was found that the static RPA 
susceptibility along c-direction in quasi-2D-tetragonal d- 
wave superconductor under strong enough magnetic field 
applied in the basal plane is singular. The instability oc- 
curs at wave-vector Q that connects points of the Fermi 
line with a finite gap. Because of the proximity of the 
compound to antiferromagnetism this can be commensu- 
rate or incommensurate close to commensuration. This 
presents evidence that the feedback of d a .2_ y 2-wave su- 
perconductivity can have drastic consequences not only 
on collective excitations in the system but also on the 
ground state properties when a external parameter tunes 
the excitation energy to zero (here realised with magnetic 
field applied in the tetragonal crystal plane). Several 
questions on this transition however remain, in particu- 
lar spectroscopic probe such as Electron Spin Resonance 
(free of low-energy incoherent peak) might help at study- 
ing closely the transition. 
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